Sensitivity coefficients to a-variation for fine-structure transitions in Carbon-fike ions 
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We calculate sensitivity coefficients to q- variation for the fine-structure transitions (1,0) and (2,1) 
within ^Pj[2s^2p^] multiplet of the Carbon-like ions C I, N ll, O ill, Na vi, Mg vil, and Si ix. These 
transitions lie in the far infrared region and are in principle observable in astrophysics for high 
redshifts z ~ 10. This makes them very promising candidates for the search for possible a- variation 
on a cosmological timescale. In such studies one of the most dangerous sources of systematic errors 
is associated with isotope shifts. We calculate isotope shifts with the help of relativistic mass shift 
operator and show that it may be significant for C I, but rapidly decreases along the isoelectronic 
sequence and becomes very small for Mg vil and Si ix. 
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A. Introduction 

Some of the extensions of the Standard Model pre- 
dict small space-time variations of such fundamental con- 
stants as the fine-structure constant a — e^jijic) and 
electron-to-proton mass ratio [i = mi, /nip. These varia- 
tions are now intensively searched for both in astrophysi- 
cal data and in laboratory experiments. For this purpose 
two, or more lines with different dependencies on the fun- 
damental constants are compared to each other at dif- 
ferent times. Laboratory experiments provide extremely 
high sensitivity to frequency shifts, while the typical 
timescale is of the order of one year. The astrophysical 
extragalactic observations allow studying time intervals 
comparable to the lifetime of the universe (~ 10^" years), 
but with much lower accuracy. Different theoretical mod- 
els predict different behavior of the constants from mono- 
tonic in time, to oscillations and sharp changes. The 
latter could take place, for example, during the transi- 
tion from matter dominated to dark energy dominated 
universe. 

We can conclude that astrophysical and laboratory 
methods are complementary and equally important. Re- 
cent progress and perspectives of the laboratory searches 
were summarized in At present all laboratory ex- 
periments give only strong upper bounds on the time- 
variation of fundamental constants. Situation in astro- 
physics is less clear. Some of the observations indicate 
non-zero variation at a few-sigma level 0, [3, 0, li while 
other results are consistent with zero variation fg, 0, Q ■ 
For more references and the discussion of the most recent 
developments see online material from the Workshop Q . 

Controversial astrophysical results may indicate some 
systematic effects. One of the possible systematic fre- 
quency shifts can arise from the different velocity dis- 
tributions of different species in molecular clouds, i.e. 



the so-called Doppler noise [l^, [ll[ . Another potentially 
dangerous systematic effect may be associated with the 
cosmological evolution of the isotope abundances, which 
can lead to the time-dependent isotope shifts of atomic 
transitions [3, • The former effect can be suppressed 
when different lines of the same specie are used to study 
possible variation of constants. The latter systematics 
is absent in microwave spectra of molecules, where lines 
of different isotopes are well resolved. Alternatively, one 
can use either atoms with single stable isotope, or specific 
combinations of atomic frequencies, which are not sensi- 
tive to isotope shifts [l^. In order to apply this latter 
method one needs to know accurate values of the isotope 
shift coefficients. 

An interesting opportunity to search for a-variation 
is associated with Ml transitions between fine-structure 
(FS) levels of ground multiplets. These lines are observ- 
able for very large cosmological redshifts, up to z ~ 10 
[fl [Tsj , and have high sensitivity to a- variation. In the 
case of ^Pj multiplet of C li ion there is only one FS 
line (1,0). As a reference one can use microwave molec- 
ular lines, which are insensitive to a-variation, but de- 
pend on the mass ratio /i. This way one can put a 
limit on the variation of the combination F = o? j [i (ll| : 
t^FjF = (O.libl.O) X 10""* at z = 6.42 (the look back time 
about 13 Gyr). The signal in this method is enhanced 
by the large and different sensitivities of compared lines. 
On the other hand, the lines of different species are used, 
so the Doppler noise can become a problem when the 
accuracy is increased. 

Many atoms and ions have triplet ground state ^Pj, 
and two Ml transitions (1,0) and (2,1) can be observed. 
The most important species of this kind for astrophysi- 
cal observations are Carbon-like and Oxygen-like ions. In 
this case one can measure the ratio of these two frequen- 
cies and compare it to the laboratory value. This way 
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the Doppler noise can be significantly suppressed. How- 
ever, in the first approximation within the LS-coupling 
scheme, the FS frequencies obey the Lande rule: 

^2,1/^1,0 = 2. (1) 

Thus, the frequency ratio does not depend on a. 

Sensitivity coefficients Q to a-variation are defined as: 

Aoj/uj = 2QAa/a. (2) 

In the approximation ([1]) these coefficients are the same 
for both FS transitions: 

Q2,i = Qi,o = 1 • (3) 

Equations ^ and ([3]) break when we take into account 
non-diagonal spin-orbit interaction. Equation ([1]) also 
breaks when Breit interaction is taken into account. If 
the latter is neglected, one can link experimentally ob- 
served violation of the Lande rule with the difference in 
sensitivity coefficients [16]: 

AQ = Q2,i-Qi.o = i f— Vl. (4) 

This simple expression predicts significant differences in 
sensitivity coefficient of FS transitions, which rapidly 
grow with nuclear charge Z. In fact, Breit interaction 
can not be neglected for the light ions with Z < 10. For 
this case Eq. ^ significantly overestimates AQ. 

For the configuration ns^np^ one can take Breit in- 
teraction into account within well known semi-empirical 
theory p7| . Corresponding results for a number of ions 
of astrophysical interest are given in [l^. However, this 
theory is essentially one-configurational and does not ac- 
count for certain correlation corrections. The latter can 
be adequately treated only within rclativistic ah initio 
calculations. 

In this paper we report ah initio calculations of the 
sensitivity coefficients Q for FS transitions in the ground 
multiplets of Carbon and Carbon-like Nitrogen, Oxygen, 
Sodium, Magnesium, and Silicon. Our results are in a 
good agreement with analytical estimates [l^ , the differ- 
ences being typically on the order of 10%. All mentioned 
elements, with exception of Sodium, have several stable 
isotopes. Inhomogeneity of isotope distribution and cos- 
mological evolution of isotope abundances can lead to 
the isotope shifts of observed lines in comparison to the 
laboratory frequencies. This can cause systematic effects 
for a-variation studies. Unfortunately, the isotope shifts 
for the FS transitions are not known. For this reason we 
calculated isotope mass shift coefficients fcms for all ions 
with several isotopes. For the ions considered here the 
volume shift is much smaller and can be neglected. 

Among the light elements C, N and O (CNO), it is 
mainly carbon, where the cosmic evolution may change 
isotopic ratios considerably. The ^^C/^^C ratio is about 
90 in the solar system, but it can be as low as 6 in the 



atmospheres of evolved massive red supergiants like Al- 
pha Ori [1^. The enhanced ^^C abundance relative to 
^^C is caused by the CNO cycle operating in massive 
stars. Therefore, one can expect that in high redshift 
objects, where only short-lived, massive stars can have 
evolved, the ^^C/^'^C ratio can be much below the solar 
value and a ^^C fraction of up to 20% is possible in the 
early universe. The question of ^^C/^^C evolution has 
also been discussed by Fenner et al. 19] for somewhat 
later phases where intermediate mass stars convert ^^C 
into ^^C in the Asymptotic Giant Brunch. Here peak 
values of -'^^C/^^C of 0.06 are predicted. For nitrogen, 
the ^''N/^^N ratio is always large so that isotopic shifts 
are probably not important. However, in the Murchison 
meteorite, the ^''N/^^N ratio varies in SiC grains between 
10^ and 4 x 10^, and the lower values of ^''N/^^N have 
been found in grains with low ^^C/^'^C [2^. For oxygen, 
the solar system value ^^O/^^O is about 2630, but this ra- 
tio can be lower in CNO processed material by a factor of 
5 (Alpha Ori, [l3|)- This means that there may be some 
cosmic evolution of isotopic ratios in N and O, however 
at a much lower level than in C, and not detectable at 
the present achievable spectral resolution. For the evo- 
lution of the isotopic ratios of Mg and Si, which have 
both several abundant isotopes, we refer to the works of 
Ashenfelter et al. [2l| and Fenner et al. [Tgj . 



B. Details of the calculations 

We used Dirac-Breit Hamiltonian and all-electron con- 
figuration interaction (CI) method. It was shown in Ref. 
|.16i] that deviation from the Lande rule was caused pri- 
marily by non-diagonal spin-orbit interaction between 
the levels ^Po,2 and ^^o, ^-D2, while the level ^Pi remained 
practically unperturbed. Therefore, it was essential that 
the theory accurately reproduced the spacings for all 5 
levels of the ground configuration. Below we present our 
results for all levels of interest and compare theoretical 
energies with the experimental data from NIST [2^ . 

Most calculations were done with Dirac-Sturmian ba- 
sis set [H, [l^!, which included orbitals up to 9s, 9p, 
8ci, 8/, and 7g. Configurational space corresponded to 
triple and partly quadruple excitations from Is- 2s- and 
2p-shells. We also did calculations with a much smaller 
configurational space which included singles for Is-shell 
and doubles for 2s- and 2p-shells, but on a somewhat 
longer basis set. In this smaller calculation we neglected 
retardation part of the Breit interaction. We found out 
that retardation corrections were very small. Expansion 
of the configurational space was more important and led 
to better agreement with the experiment for the transi- 
tion frequencies. On the other hand, the difference for 
the sensitivity coefficients was not dramatic (see below). 
Correlation effects were most important for neutral Car- 
bon and less important for the Carbon- like ions. 

Comparison of the results for frequencies and Q-factors 
from the calculations with two configurational spaces de- 
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TABLE L Energies in cm ^ and Q-factors for the levels of the configuration 2s^2p^ in respect to the ground state ^Pq. 
Experimental frequencies are taken from Ref. [2^ . 



Ion 


Exp Gr . 


■in 
ri 

Theor. 


Q 


Expcr. 


Theor. 


Q 


Exp Gr . 


Theor. 


Q 


Expcr. 


Of) 

Theor. 


Q 


C I 


16.42 


16.4 


0.9998 


43.41 


43.3 


0.9948 


10192.6 


10435 


0.0023 


21648.0 


21879 


0.0014 


N II 


48.7 


48.3 


1.0086 


130.8 


129.3 


1.0004 


15316.2 


15605 


0.0052 


32688.8 


33036 


0.0031 


O III 


113.2 


112.8 


1.0197 


306.2 


304.9 


1.0040 


20273.3 


20535 


0.0099 


43185.7 


43540 


0.0058 


Na VI 


698 


693 


1.0671 


1859 


1850 


1.0132 


35498 


35719 


0.0375 


74414 


74795 


0.0215 


Mg VII 


1107 


1108 


1.0891 


2924 


2920 


1.0172 


40948 


41192 


0.0534 


85153 


85571 


0.0306 


Si IX 


2545 


2532 


1.1403 


6414 


6394 


1.0254 


52926 


53185 


0.0990 


107799 


108253 


0.0573 



scribed above confirmed that the smaller space was suf- 
ficiently saturated. That allowed us to use it for the 
calculations of the isotope shifts and significantly reduce 
the computational costs. 

In order to determine Q-factors we calculated atomic 
energy levels for three or five values of a in the vicinity 
of its physical value and performed numerical differenti- 
ation. Three point differentiation was used in most pre- 
vious calculations [l^, [2^. Here we were interested in 
the small differences in Q-factors for the FS levels, so we 
compared three-point differentiation with more expensive 
five-point differentiation. The difference appeared to be 
rather small. 

Calculations of the mass shift coefficients are similar to 
calculations of Q-factors [13, HI]. The Hamiltonian iJms, 
which describes mass shift, is added to atomic Hamilto- 
nian with the coefficient A: il\ = H + Affms- After solv- 
ing eigenvalue problem for H\ one finds the mass shift 
coefficient by numerical differentiation, fcms = dE\/dX. 

For the mass shift we use relativistic theory developed 
m For optical transitions in light atoms 

the relativistic corrections to mass shift are usually small 
[S^ . However, here we are interested in the infrared FS 
transitions. As long as the fine-structure disappears in 
the non-rclativistic limit, the non- relativistic theory of 
the isotope shift is not applicable here. 

In the non-relativistic approximation the total mass 
shift is a sum of the normal mass shift (NMS) and the 
specific mass shift (SMS). The NMS contribution comes 
from the substitution of the electron mass with the re- 
duced mass, and respective coefficient is simply propor- 
tional to the transition frequency: k™^^ = a;/1823 (a.u.). 
The SMS contribution is described by the two-electron 
operator J2i^k Pi ' Pk has to be calculated numeri- 
cally. In the relativistic theory the NMS and SMS con- 
tributions correspond to one-electron (i = k) and two- 
electron (i ^ k) parts of the relativistic mass shift oper- 
ator: 

Hms = X! ~ + ■ '^)^]) 

i,k ^ ' i 

'{p~^ \oL + {OL ■ r)r\\^ , (5) 



which is correct to the second order in aZ. Here NMS 
contributions have to be calculated numerically on the 
same footing as SMS ^29,, i3^, [Ml, [H [H . 

C. Numerical results and discussion 

Results of our calculations of the energies and sensi- 
tivity coefficients Q for the levels of the ground config- 
uration ls^2s^2p^ are presented in Table [H All of them 
are given in respect to the ground state ■^Pq- Calculations 
are done for the "large" configurational space described 
in Sec. \B\ This space includes ~ 1.2 x 10^ relativistic 
configurations. 

We see that for the levels of the ground multiplet the 
Q-factors are close to unity, as expected. For two other 
levels of the ground configuration the Q-factors are small. 
Deviation from unity in the former case and deviation 
from zero in the latter are of the same order of magnitude 
and rapidly grow with nuclear charge Z. Here we are 
mostly interested in the differences of the Q-factors for 
the FS transitions (1-0) and (2-1). These differences are 
given in Table HIl 

From Table |TT] one can see that our numerical results 
are in goocl_ agreement with semi-empirical estimates 
from Ref. [l(|. Calculations with "large" and "small" 
CI described above are very close to each other. Except 
for the neutral Carbon, the difference between two nu- 
merical calculations is much smaller than their difference 
from semi-empirical values. Note, that for the large CI 
we used five-point differentiation method and three-point 
method for the small CI. Therefore, we can conclude that 
small CI is already sufficiently saturated and that three- 
point differentiation is also sufficiently accurate. Because 
of that, in the following calculations of the mass shifts we 
use this significantly cheaper method (small CI includes 
~ 10000 relativistic configurations). 

Our results for the mass shifts are summarized in Ta- 
bles mil and IIVI Table IIIII gives coefficients fcnms and 
^sms, while Table [TVl presents full mass shifts for given 
isotope pairs. We see that for the FS transitions from 
the ground state the NMS coefficients calculated with 
the help of relativistic operator ([5]) differ very strongly 
from the non-relativistic values. Interestingly, the differ- 
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TABLE 11: The differences of the sensitivity coefficients AQ of the FS emission hnes within the ground multiplet '^Pj for the 
most abundant Carbon-hke ions. Ab initio results of this work are compared with semi-empirical results [l6l ] based on the 
theory [13]. In addition to the large CI we also made a much smaller CI as described in the text. 

Ion Transition a (1,0) Transition b (2,1) uJb/i^a AQ = Qb — Qa 





K (^ni) 


LUa (cm ^) 


Af, (^m) 


LUb (cm ^) 




Ref.[16] 


Large CI 


Small CI 


C I 


609.1 


16.40 


370.4 


27.00 


1.646 


-0.008 


-0.0081 


-0.0090 


N II 


205.3 


48.70 


121.8 


82.10 


1.686 


-0.016 


-0.0132 


-0.0134 


III 


88.4 


113.18 


51.8 


193.00 


1.705 


-0.027 


-0.0250 


-0.0247 


Na VI 


14.3 


698 


8.6 


1161 


1.663 


-0.091 


-0.0861 


-0.0843 


Mg VII 


9.0 


1107 


5.5 


1817 


1.641 


-0.12 


-0.116 


-0.114 


Si IX 


3.9 


2545.0 


2.6 


3869 


1.520 


-0.21 


-0.190 


-0.188 



TABLE III: Coefficients fcnms and fcsms for Carbon-like ions. 
Mass shift is calculated in respect to the ground state 
All numbers are in (GHz x amu). For NMS we also give non- 



relativistic value fc" 



= A£e: 



r(a.u.) X 3609.6. 







'P. 


'P2 






C I 


^nms 


-0.09 


-0.42 


178.8 


369.7 




"'nms 


0.27 


0.71 


167.6 


356.0 




^sms 


0.51 


1.75 


-155.2 


-151.2 


N II 


h 


-0.31 


-1.26 


263.5 


553.2 




"'nms 


0.80 


2.15 


251.9 


537.6 




h 


1.35 


4.42 


-179.4 


-167.4 


III 


^nms 


-0.09 


-1.43 


343.3 


725.2 




"'nms 


1.86 


5.04 


333.4 


710.2 




^sms 


2.56 


8.79 


-197.0 


-165.6 


Mg VII 


^nms 


7.33 


14.51 


661.9 


1398.3 






18.21 


48.09 


673.4 


1400.4 




^sms 


21.58 


62.86 


-221.5 


131.1 


Si IX 


^nms 


19.83 


43.96 


834.0 


1748.3 




'^nms 


41.86 


105.49 


870.4 


1772.9 




^sms 


45.34 


124.24 


-197.6 


468.3 



TABLE IV: Frequency shifts 5v — — for FS transitions 
(J' — J) in isotopes A' and A of C-like ions. The last two 
columns give respective velocity shifts SV — —Su/i/ x c, where 
c is the speed of light. 



Siy (GHz) 



5V (km/s) 



Ion 


A' 


A 


(0-1) 


(1-2) 


(0-1) 


(1-2) 


C I 


13 


12 


-0.00278 


-0.00602 


1.691 


2.228 


N II 


15 


14 


-0.00498 


-0.01006 


1.024 


1.225 


III 


17 


16 


-0.00909 


-0.01798 


0.803 


0.932 




18 


16 


-0.01717 


-0.03396 


1.517 


1.760 


Mg VII 


25 


24 


-0.04819 


-0.08076 


0.435 


0.444 




26 


24 


-0.09267 


-0.15531 


0.837 


0.855 


Si IX 


29 


28 


-0.08026 


-0.12689 


0.315 


0.328 




30 


28 


-0.15517 


-0.24532 


0.610 


0.634 



ence is bigger for the light ions, where even the sign of 
the non-relativistic approximation is wrong. Relativistic 
NMS coefficient becomes positive between Z = 8 (Oxy- 
gen) and Z = 9 (Fluorine), but remains significantly 



smaller than non-relativistic predictions. Note that for 
two optical transitions between levels of different multi- 
plets the relativistic and non-relativistic values are very 
close. This result confirms that non-relativistic theory of 
the mass shift works nicely for optical transitions, but is 
not applicable to FS transitions. 

Table Hill shows that SMS coefficients also significantly 
change along the isoelectronic sequence. The NMS coef- 
ficients are mostly sensitive to relativistic effects, which 
depend on the parameter aZ. The SMS coefficients are 
more sensitive to correlations. The latter are governed 
by the parameter l/Z and decrease along the sequence. 
It is clearly seen from Table IIVI that full mass shift for 
the FS transitions behaves more smoothly, than NMS 
and SMS parts. This may mean that in the relativistic 
theory there is no good reason to separate mass shift into 
NMS and SMS parts. The last two columns of Table HVl 
give apparent velocity (Dopplcr) shifts in astrophysical 
observations, which correspond to the frequency shifts 
from two previous columns. We see that velocity shifts 
for two FS transitions significantly differ for light ions, 
but become almost equal for heavier ions. For Mg vii 
and Si ix the isotope shifts practically cancel out from 
the frequency ratio W2,i/wi,o- Consequently, for these 
ions the isotope shift does not lead to noticeable system- 
atic effect in a- variations studies. 



D. Conclusions 

FS transitions can be observed in far infrared wave- 
band for very distant astrophysical objects with redshifts 
z ~ 10. This makes them very promising candidates 
for the study of possible a-variation on the cosmologi- 
cal timescale. We performed ab initio calculation of the 
sensitivity coefficients Q to a-variation for the FS tran- 
sitions within ground multiplet '^Pj of Carbon-like light 
ions. These calculations confirmed that the differences in 
sensitivity for FS transitions A Q is of the same order as 
AQ for optical transitions in the same ion. In both cases 
AQ rapidly grows with Z. In optical waveband the most 
dangerous systematic effect is associated with the isotope 
shifts. We calculated mass shifts for FS transitions and 
found that isotope shift is rather large for light ions C i. 
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N II, and O ill, but practically cancels out for Mg Vii and 
Si IX. At the same time these heavier ions, together with 
the Na vi ion, which has only one stable isotope, have 
higher sensitivity to a- variation. 

It was pointed out in Sec. El that in the solar sys- 
tem light elements C, N, and O have only one dominant 
isotope. However, in the early universe the abundances 
of ^^C and ^^C could be significantly different. At the 
same time, the isotope shift for the FS transitions in 
C I is larger than for other C-like ions, while sensitiv- 
ity to a-variation is lower. We conclude that the method 
suggested in [l6| may not work for C i because of the 
small sensitivity and large systematic effects associated 
with isotope shifts. For C-like ions the sensitivity to a- 
variation grows with the nuclear charge Z, while isotope 
shifts rapidly decrease. Moreover, for heavier ions the 



isotope shifts almost cancel out for the ratio of the FS 
frequencies. Another dangerous systematic effect from 
the Doppler noise is significantly suppressed for the pairs 
of lines of the same species. Therefore, observations of 
the FS transitions in heavier C-like ions can be used as 
a sensitive tool for the search of a-variation at large red- 
shifts. 
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